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RESEARCH  GOALS: 


During  Apri 1-May- June  1987  and  June-July-August  1988  the 
Coastal  Transition  Zone  (CTZ)  Experiments  were  conducted  off  the 
California  Coast  near  Pt.  Arena  by  the  Office  of  Naval  Research 
(ONR)  .  One  of  the  basic  goals  of  these  experiments  was  :  "To 
understand  the  role  of  the  atmosphere  in  driving  coastal  ocean 
currents  and  to  learn  of  possible  ocean-atmosphere  feedbacks." 
This  component  of  the  CTZ  program  was  the  focus  of  the 
meteorological  conditions  during  these  experiments  including  the 
data  collection,  reduction,  and  interpretation.  This  grant  covered 
primarily  the  1988  study. 

OBJECTIVES : 


This  principal  investigator  acted  as  a  focal  person  for  the 
meteorology  of  the  CTZ  experiments  providing  meteorological 
expertise  to  the  other  investigators.  Our  basic  objectives  were: 


1.  Documentation  of  the  meteorological  conditions  during 
each  field  study. 


•2. 


Data  reduction  plus  analysis-synthesis  of  the 
meteorological  conditions  for  each  field  study. 


Yrm  quality  inspected  3 
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ACCOMPLISHMENTS: 


This  P.I.  was  not  in  the  field  during  the  CTZ  experiments,  but 
provided  weather  outlooks  and  wind  data  to  the  P.I.'s  via  OMNET 
during  the  experiments.  Our  accomplishments  are  listed  below. 

DOCUMENTATION:  We  acquired  extensive  holdings  of  the 
meteorological  data  for  the  field  phase  and  these  data  have  been 
made  available  to  other  P.I.'s  as  needed.  Once  per  day  (OOGMT) 
surface  weather  maps  are  on  file  for  the  period  of  1  June-31  August 
1988.  Loops  of  satellite  cloud  images  for  the  CTZ  area  for  the 
period  of  1  June-23  August  1988  via  GOES-West  are  readily  available 
on  video  tapes.  The  National  Weather  Service  (NWS)  daily  show 
known  as  A.M.  Weather  is  also  available  for  the  period  of  1  June-26 
August.  Twice  weekly  maps  of  SST,  from  a  composite  of  data 
sources,  prepared  by  NWS  at  their  San  Francisco  Forecast  Office  are 
in  our  archives  for  the  period  of  7  June-30  August  1988  in  hard 
copy  map  format.  Also  hard  copies  of  hourly  weather  reports  for 
the  West  Coast  are  in  our  files  for  the  period  28  June-12  August 
1988. 

Our  documentation  is  very  extensive  for  the  National  Data  Buoy 
Center's  (NDBC)  coastal  buoys  and  C-MAN  stations.  We  have  hard 
copies  of  the  hourly  data  for  the  period  1  March-31  August  for 
Buoys  46013  (Bodega  Bay)  and  46014  (Pt.  Arena)  and  C-MAN  PTAC1  (Pt. 
Arena) .  We  also  have  these  and  all  the  West  Coast  Buoy  and  C-MAN 
station  data  on  magnetic  tape  for  the  period  1  March-31  August 
1988. 
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These  buoy  and  C-MAN  stations  are  shown  in  Fig.  1  plus  NDBC-05 
(Washington) ,  28  (Cape  San  Martin) ,  40  (Cape  Foulweather) ,  42 
(Monterey)  and  C-MAN  CAR03  (Cape  Arago)  and  NWP03  (Newport) . 

In  short,  we  were  quite  successful  in  obtaining  the  standard 
data  via  the  standard  station  spacing  to  document  the  broad 
meteorological  conditions  during  the  field  phase  as  well  as  before 
and  afterwards.  Many  of  these  data  were  shared  with  the  CTZ  P.I.'s 
during  and  after  the  experiment. 

DATA  PLOTS:  The  NDBC  buoy  and  C-MAN  data  have  been  prepared 
into  monthly  time  plots  to  help  show  the  response  of  Sea  Surface 
Temperature  (SST)  to  the  local  wind  stress.  These  time  plots  use 
hourly  data  but  after  it  has  been  smoothed  using  a  25  point  (hour) 
running  unweighted  mean.  We  have  prepared  these  monthly  time  plots 
for  all  of  the  buoys  and  C-MAN  stations  cited  above.  Plots  were 
prepared  for  the  sea  level  pressure,  air  temperature,  sea  surface 
temperature  (SST) ,  wind  direction  and  speed,  u  (E-W)  and  v  (N-S) 
components  of  the  wind,  and  wind  stress  (magnitude  and  direction) 
for  each  of  these  stations.  We  include  a  composite  of  stress 
(magnitude  and  direction)  and  SST  for  NDBC  Buoy  13  (Bodega  Bay)  for 
June  July  August  1988  (Figure  2,3,4)  to  illustrate  the  data  format 
made  available  to  the  P.I.'s  during  the  synthesis  stage  of  this 
work. 

These  plots  show  clear  examples  of  SST  response  to  increases 
and  decreases  in  the  wind  stress.  This  is  especially  true  in  the 
June  and  July  plots,  but  it  is  also  evident  in  the  August  plots. 
Such  graphs  were  very  valuable  at  the  CTZ  workshop  in  Corvallis 
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during  25-27  July.  The  "Filament  Tracking  Group"  focused  on  the 
period  4-11  July  when  there  were  three  Research  Vessels  (WECOMA, 
Pt.  Sur,  and  Washington)  in  the  area.  The  "Change  in  Filament 
Group"  focused  on  a  period  in  late  July. 

MANUSCRIPTS :  Many  of  the  research  results  of  the  CTZ 
Experiments  in  1987  and  1988  appear  in  a  special  section  of  the 
August  15,  1991  issue  of  the  Journal  of  Geophysical  Research.  This 
P.I.  worked  with  many  of  the  other  P.I.'s  particularly  in  the  two 
groups  mentioned  above.  The  wind  data  and  general  meteorological 
conditions  are  woven  throughout  the  numerous  papers.  However,  this 
P.I.  was  a  coauthor  only  on  the  paper  by  Brink  et  al.  A  reprint  of 
that  paper  is  attached. 

MEETINGS:  This  P.I.  attended  the  CTZ  Workshops  in  September 
1987  in  Corvallis,  OR,  in  November  1987  in  LaJolla,  CA,  and  again 
in  July  1989  in  Corvallis,  OR.  This  latter  meeting  was  during  this 
grant  period  and  was  the  meeting  at  which  the  groups  cited  above 
were  organized  and  the  framework  of  the  manuscripts  took  place. 
The  P.I.  also  attended  the  American  Meteorological  Society  meeting 
in  New  Orleans  in  January  1991  with  visits  to  the  ONR  Laboratories 
in  Bay  St.  Louis,  MS. 

BUDGET:  The  main  employees  under  this  grant  were  the  P.I. 
(for  limited  support  during  the  Summers  of  1989,  1990,  and  1991), 
Robert  Joyce  (graduate  student-Summer  1989) ,  and  Pejaver  Rao  (under 
graduate  student -Fall  1991)  .  Expense  items  covered  minor  supplies, 
phone  services  including  OMNET,  postage,  xeroxing  plus  computer 
supplies . 


5 


The  total  funds  available  were  $39,  995  over  the  full  period  of  this 
grant . 

DISTRIBUTION  LIST: 

This  final  technical  report  is  supplied  as  follows: 


ADDRESSEES  NUMBER  OF  COPIES 

Scientific  Officer  Code:  1121CS  3 

Thomas  Kinder 

Office  of  Naval  Research 

800  North  Quincy  Street 

Arlington,  Virginia  22217-5000 

Administrative  Grants  Officer  1 

Office  of  Naval  Research 

Resident  Representative  N66020 

Administrative  Contracting  Officer 

Georgia  Institute  of  Technology 

206  O'Keefe  Building 

Atlanta,  Georgia  30332-0490 

Director,  Naval  Research  Laboratory  1 

Attn:  Code  2627 

Washington,  D.C.  20357 

Defense  Technical  Information  Center  2 

Building  5,  Cameron  Station 
Alexandria,  Virginia  22314 

Robert  M.  Johnson,  Vice  President  1 

Sponsored  Research 
Florida  State  University 
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Figure  1  -  Meteorological  Stations  in  the  CTZ  Area 
Data  plots  were  prepared  for  the  underlined  buoy 
and  C-MAN  stations 


SEA  SFC  TE HP  (Cl  WIND  STRESS  (DYNES/CM»CM)  STRESS  -  DIRECTION  IOEG1 


BUOY  13  BOOCCA  MV  JUNE  IBM 


SMOOTHED 


TIME  CGHT-OAVSI 


Figure  2  -  Time  Plots  for  NDBC  Buoy  13  Bodega  Bay  -  June  1988 
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Figure  4  -  Time  Plots  for  NDBC  Buoy  13  Bodega  Bay  -  August  1988 
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Statistical  Properties  of  Near-Surface  Flow  in  the  California 

Coastal  Transition  Zone 

K.  H.  Brink,1  R.  C.  Beardsley,1  P.  P.  Njjler,3  M.  Abbott,3  A.  Huyer 3  S.  Ramp,4  T.  Stanton  4 

and  D.  Stuart* 


During  the  iiunmen  of  1967  and  1988,  77  near-surface  satellite-tracked  drifter*  were  deployed 
in  or  near  cold  filament*  near  Point  Arena,  California  (39*N),  and  tracked  for  up  to  6  month, 
a*  part  of  the  Coastal  Transition  Zone  (CTZ)  program.  The  drifter*  had  large  drogues  centered 
at  IS  m,  and  the  resulting  drifter  trajectory  data  set  has  been  analysed  in  terms  of  its  Eulerian 
and  Lagrangian  statistics.  The  CTZ  drifter  result*  show  that  the  California  Current  can  be 
characterised  in  summer  and  fall  as  a  meandering  coherent  jet  which  on  average  Sows  southward 
to  at  least  30* N,  the  southern  end  of  the  study  domain.  FVom  39*N  south  to  about  33'N,  the 
typical  core  velocities  are  of  0(50  cmi"1)  and  the  current  meanders  have  alongshore  wavelength* 
of  0(300  km)  and  onshore-ofbhore  amplitude  of  0(100-300  km).  The  lateral  movement  of  this 
jet  leads  to  large  eddy  kinetic  energies  and  large  eddy  diffusivities,  especially  north  of  36*N.  The 
initial  onshore-offshore  component  of  diffusivity  is  always  greater  than  the  alongshore  component 
in  the  study  domain,  but  at  the  southern  end,  the  eddy  diffusivity  is  more  isotropic,  with  scalar 
single  particle  diffusivity  (KXI  +  K**)  of  0(8  X  107  cm*  s-1).  The  eddy  diffusivity  increases 
with  increasing  eddy  energy.  Finally,  a  simple  volume  budget  for  the  1988  filament  observed 
near  37*  N  off  Point  Arena  suggests  that  lubduction  can  occur  in  a  filament  at  an  average  rate 
of  0(10  m  d-1)  some  300  km  offshore,  thus  allowing  the  cold  water  initially  in  the  filament 
core  to  sink  below  the  warmer  ambient  water  by  the  time  the  surface  velocity  core  has 
turned  back  onshore.  This  process  explains  why  satellite  temperature  and  color  imagery  tend 


to  “see”  only  flow  proceeding  offshore. 


1.  Introduction 

One  central  goal  of  the  Coastal  Transition  Zone  (CTZ) 
program  has  been  to  understand  the  filaments  of  cold  water 
which  project  offshore  from  the  coast  of  northern  California 
during  the  summer  upwelling  season.  To  help  achieve  this 
goal,  intensive  physical  oceanographic  measurements  were 
made  in  the  upper  ocean  off  California  in  the  summers  of 
1987  and  1988  to  provide  new  and  more  complete  descrip¬ 
tions  of  the  water  property  and  velocity  structure  of  these 
cold  filaments  and  their  surrounding  environment.  Based  in 
part  on  the  successful  early  exploration  into  the  nature  of  fil¬ 
aments  by  Davis  [1985a]  using  radio-tracked  surface  drifters, 
the  CTZ  field  program  included  a  large  satellite- tracked 
drifter  program  to  observe  and  characterize  the  near-surface 
flow  and  its  dynamical  properties.  Since  the  filaments  were 
known  to  be  surface-intensified  features,  the  CTZ  drifters 
were  all  drogued  at  a  common  depth  of  15  m,  which  fre¬ 
quently  placed  the  drogue  within  the  surface  mixed  layer. 

The  objective  of  this  study  is  to  use  the  CTZ  drifter  mea¬ 
surements  to  characterize  statistically  the  kinematics  anu 
dynamics  of  the  near-surface  (15  m)  flow.  Since  drifters 
tend  to  spread  out  over  a  broad  area  with  time,  making 
full  use  of  the  drifter  data  base  often  implies  looking  on  a 
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much  broader  spatial  scale  than  that  encompassed  by  the 
rest  of  the  CTZ  measurements,  which  were  confined  to  the 
area  north  of  roughly  37*N  (Figure  1).  The  drifter  data 
used  here  cover  latitudes  from  about  30*  to  39*N,  that  is 
to  say,  most  of  the  coast  of  California.  Such  broad  scale 
descriptions  will  sometimes  tend  to  deemphasize  specific  fil¬ 
amentary  features,  although  these  are  of  considerable  inter¬ 
est.  A  description  of  the  individual  CTZ  drifter  trajectories, 
including  the  tendency  for  the  near-surface  (15  m)  flow  to 
proceed  offshore  along  the  axis  of  a  cold  filaments,  is  pre¬ 
sented  by  M.  Swenson  et  al.  (Drifter  observations  of  the 
dynamical  and  thermodynamical  structures  in  a  cold  fila¬ 
ment  off  Point  Arena,  California,  in  July  1988,  submitted  to 
Journal  of  Geophysical  Research,  1990,  hereinafter  referred 
to  as  Swenson  et  al.,  1990). 

The  following  results  deal  with  a  number  of  specific  top¬ 
ics,  all  addressed  statistically.  First,  the  “Eulerian*  mean 
flow  and  its  accompanying  variability  will  be  treated.  Next, 
the  dispersion  properties  measured  using  the  drifters  will  be 
discussed.  Both  of  these  first  two  topics  tend  to  deal  with 
the  broader  spatial  scales.  Next,  the  momentum  balance  for 
near-surface  flow  will  be  treated  in  the  CTZ  intensive  area. 
Finally,  the  relation  of  drifter  kinematics  to  frontal  locations 
in  the  Point  Arena  vicinity  will  be  considered. 

2.  The  Observational  Base 

The  primary  measurement  system  is  the  Tristar  Mark  II 
drifter  [Afttfer  et  at.,  1987],  These  instruments  have  large 
(6  m)  radar  reflector-shaped  drogues  located  at  a  central 
depth  of  15  m.  The  drogue  and  surface  float  design  has 
been  calibrated  by  in  situ  current  measurements  and  follows 
horizontal  currents  to  within  about  1  cm/s  under  conditions 
representative  of  those  encountered  here.  The  drifters  are 
tracked  using  the  Argos  satellite  system.  Typically,  eight 
fixes  per  day  were  obtained  during  the  CTZ  program,  and 
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Fig.  1.  Locator  map.  The  open  circle  denotes  the  location  of 
NDBC  buoy  46013.  The  solid  circles  are  locations  of  hydrographic 
stations  for  a  typical  1988  survey,  defining  the  CTZ  intensive  re¬ 
gion.  Boxes  represent  the  areas  within  which  diffiisivities  and 
Lagrangian  time  scales  were  computed.  Crosses  represent  the  lo¬ 
cations  of  the  Stabcno  and  Smith  [1987]  moorings  used  to  estimate 
Eulerian  time  scales. 


the  accuracy  per  fix  is  generally  about  300  m.  After  initial 
editing  of  the  drifter  positions,  cubic  splines  (each  4  days  in 
length)  are  used  to  smooth  the  tracks  and  remove  inertial 
and  tidal  variance.  The  splines  effectively  represent  a  low- 
pass  filter  with  a  half  power  point  of  about  0.45  cpd. 

The  drifters  were  deployed  during  the  late  spring  or  sum¬ 
mer  of  1987  and  1988  in  small  coherent  clusters.  Most  of  tne 
25  deployments  in  1987  were  in  May  and  June,  while  most 
of  the  52  in  1988  were  in  July  or  early  August.  In  1987 
the  drifters  provided  good  data  for  typically  25  days,  while 
in  1988,  good  data  were  typically  available  for  6  months  or 
more  following  launch.  The  drifters  were  usually  launched 
in  clusters  of  3-17,  with  separations  of  5-20  km  from  near¬ 
est  neighbors.  Smaller  separations  were  used  in  larger  de¬ 
ployments.  The  deployment  strategy  was  intended  to  al¬ 
low  mapping  of  velocity  and  relative  vortidty  within  the 
Point  Arena  filament  (Swenson  et  al.,  1990).  The  total  set 
of  drifter  tracks  (Figure  2)  gives  a  measure  of  the  amount  of 
data  available. 

In  addition  to  the  drifter  trajectories,  several  other 
sources  of  data  were  employed.  Conductivity- temperature- 
depth  (CTD)  surveys  were  made  repeatedly  on  a  regular  grid 
(Figure  1)  in  the  CTZ  intensive  area.  Details  on  the  hydro- 
graphic  observations  are  given  by  Fleischbein  et  al.  (1989] 
and  Huyer  et  al.  [this  issue].  Winds  were  measured  at  10  m 
elevation  at  NOAA  Data  Buoy  Center  (NDBC)  environmen¬ 
tal  buoy  46013,  located  at  38*13'N,  123#18'W.  The  wind 
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Fig.  2.  Complete  (a)  1987  and  (1)  1988  CTZ  drifter  tradu. 
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velocity  time  series  were  converted  into  wind  stresses  using 
the  method  of  Large  and  Pond  [1981],  and  then  spline- fit  in 
the  same  way  as  the  drifter  positions  to  provide  a  low-pass- 
filtered  wind  stress  time  series.  Finally,  advanced  very  high 
resolution  radiometer  (AVHRR)  images  of  surface  tempera¬ 
ture  were  used.  These  were  calibrated  and  registered  using 
the  techniques  described  by  Kelly  and  Davit  [1986],  except 
that  no  view  angle  correction  was  made. 

3.  Eulerian  Statistics 

The  general  properties  of  the  near-surface  flow  will  first  be 
described  in  terms  of  “Eulerian”  flow  patterns.  We  use  the 
term  “Eulerian”  cautiously  because  of  the  bias  introduced 
by  the  tendency  for  drifters  to  avoid  divergences  and  to  enter 
an  area  only  when  the  flow  is  in  the  right  sense.  In  order 
to  make  these  estimates,  1-day  averages  of  smoothed  drifter 
velocities  were  sorted  into  0.5*  by  0.5*  latitude-longitude 
bins.  Bins  containing  less  than  10  drifter  days  of  data  were 
deleted  from  further  consideration.  For  each  bin,  the  vector 
mean  flow  and  the  standard  deviations  along  principal  axes 
directions  were  computed.  The  results  are  shown  in  Figure  3 
for  1987  and  1988  separately. 

In  the  bins  where  there  was  a  large  data  base  (greater 
than  50  drifter  days),  standard  deviations  of  currents  were 
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Fig.  3.  Eulerian  statistics  from  drifters.  Results  are  only  shown 
for  boxes  with  more  than  10  drifter  days  of  data.  Dots  represent 
box  centers,  (s)  Mean  velocity  for  1987.  (1)  Standard  deviations 
in  principal  axis  system  for  1987.  The  half-length  of  each  cross 
represents  1  standard  deviation,  (c)  Number  o t  drifter  days  per 
box  for  1987.  (4)  Mean  velocity  for  1988.  (e)  Standard  deviations 
in  principal  axis  system  for  1968.  (/)  Number  of  drifter  days  per 
box  for  1988. 
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typically  about  20-30  cm  s-1.  The  Eulerian  independence 
time  scale  was  calculated  from  moored  current  meter  records 
in  the  region  (Figure  1)  at  depths  of  145-200  m  [Staheno 
and  Smith,  1987].  Although  these  instruments  are  far  from 
the  surface,  they  are  the  shallowest  long  records  available. 
Although  currents  in  the  region  are  definitely  baroclinic,  re¬ 
versals  of  direction  in  the  upper  200  m  generally  occur  only 
when  near-surface  currents  are  weak  («£10  cm  s-1),  [e.g., 
Kosro  et  at.,  this  issue],  so  we  believe  that  the  145-200  m 
Eulerian  current  independent  time  scales  are  at  least  rep¬ 
resentative  of  the  near-surface  flow.  The  average  value  ob¬ 
tained  from  these  records  was  12  days  (obtained  by  integrat¬ 
ing  the  autocorrelation  over  only  positive  lags  [e.g.,  Kundu 
and  Allen,  1976].  The  drifter  bins  with  the  best  data  cover¬ 
age  span  about  100  calendar  days,  or  represent  about  4  de¬ 
grees  of  freedom  using  the  current  meter  time  scale.  Thus 
the  90%  error  bounds  on  velocity  in  the  most  data-rich  bins 
are  typically  about  ±26  cm  s-1  (about  1.5  times  the  length 
of  an  arrowhead  in  Figure  3).  The  most  data-rich  bins  (Fig¬ 
ures  3c  and  3 f)  tend  to  be  those  along  the  core  of  the  mean 
jets  visible  in  Figure  3.  Note  that  in  a  sense,  this  is  a  con¬ 
servative  error  estimate,  since  it  is  based  on  calendar  days 
of  occupation  rather  than  drifter  days  in  a  box. 

The  means  from  both  years  (Figure  3)  and  from  the  two 
years  combined  (Figure  4)  show  a  meandering  core  current 
which,  on  the  whole,  proceeds  southward.  Note  that,  de¬ 
spite  similar  release  points,  the  1987  and  1988  data  sets 
have  very  different  spatial  coverages,  with  the  1987  drifters 
being  more  concentrated  inshore.  In  both  years,  the  aver¬ 
age  core  velocities  are  about  50  cm  s-1 ,  and  the  meandering 
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Fig.  4.  A*  in  Figure  3,  but  based  on  the  merged  1987  and  1988 
data  **ts.  (a)  Mean  velocity.  (1)  Standard  deviations. 


patterns  have  alongshore  wavelengths  of  0(300  km).  The 
locations  of  the  cores  differ  south  of  about  36.5*N  between 
the  two  years,  and  the  differences  are  reflected  in  data  dis¬ 
tributions  which  barely  overlap  south  of  35*N.  For  this  rea¬ 
son,  only  in  the  area  north  of  about  36*  N  can  the  1987  and 
1988  data  be  compared  meaningfully.  In  1987,  the  offshore- 
directed  branch  near  37*N  does  not  proceed  as  far  offshore 
as  the  similar  feature  near  37.5*N  in  1988.  This  is  qualita¬ 
tively  consistent  with  AVHRR  images,  which  suggest  that 
the  filaments  were  shorter  in  1987  than  in  1988.  Both  years 
show  significant  eastward  velocity  at  about  36*  N,  suggest¬ 
ing  that  at  least  this  part  of  the  meander  field  may  be  geo¬ 
graphically  stable.  It  is  possible  that  the  differences  between 
the  1987  and  1988  mean  current  fields  are  related  primar¬ 
ily  to  seasonal  variability  associated  with  the  slightly  earlier 
1987  launches.  Considerable  caution,  of  course,  must  be 
used  in  interpreting  the  mean  flow  results  because  the  Eu- 
lerian  sampling  grid  is  dependent  upon  the  flow  itself  and 
the  specific,  tightly  clustered  deployment  scheme  used  in 
the  CTZ  program.  A  further  complicating  factor  is  that, 
while  the  drifters  were  deployed  in  May-August,  most  of 
the  fixes  in  the  southern  region  were  obtained  in  Novem¬ 
ber  and  December.  Thus  there  is  a  possibility  that  sea¬ 
sonal  variations  in  the  California  Current  system  can  be 
aliased  into  apparent  spatial  variations  of  the  mean  flow. 
Climatological  0/500  dbar  flow  patterns  in  this  area  are  ori¬ 
ented  northwest-southeast  [Hickey,  1979],  In  the  summer, 
flow  is  southeastward  everywhere,  while  in  the  winter,  flow 
within  100-200  km  of  shore  reverses  to  flow  northwestward. 
In  both  seasons,  maximum  mean  alongshore  velocities  are 
about  15  cm  s-1,  and  5  cm  s_1  is  typical.  Thus  a  seasonal 
aliasing  of  our  mean  flow  seems  possible,  since  the  clima¬ 
tological  “California  Current”  as  defined  by  the  0/500  dbar 
dynamic  topography  [Hickey,  1979]  is  indeed  weaker  in  the 
winter,  although  the  combined  mean  drifter  flow  (Figure  4) 
is  both  stronger  and  directed  more  southward  than  the  cli¬ 
matological  version  [Hickey,  1979]  over  the  entire  domain. 

The  standard  deviations  of  velocity  present  a  rich  pattern. 
Perhaps  the  most  obvious  aspect  is  a  tendency  for  variance 
to  decrease  from  north  to  south  (Figure  4).  This  could  re¬ 
flect  a  true  north-south  gradient  of  the  eddy  kinetic  energy, 
seasonal  variations  in  velocity  variance  or  it  could  represent 
a  tendency  for  more  evenly  spread  velocity  sampling  farther 
to  the  south,  since  most  of  the  drifters  were  deployed  in  or 
near  energetic  structures  near  38*-39*N.  A  further  tendency 
is  for  the  major  axis  of  variability  to  line  up  loosely  with 
the  mean  flow  in  the  jet  core  regions,  most  notably  in  the 
northern  part  of  the  domain.  South  of  about  33*N,  where 
the  mean  jet  is  no  longer  very  pronounced,  the  principal  axis 
orientations  appear  to  be  fairly  random,  although  there  is  a 
tendency  for  the  axes  to  tilt  toward  the  jet  core,  suggesting 
a  convergence  of  southward  momentum.  The  tendency  for 
current  variability  to  decrease  toward  the  south  could  repre¬ 
sent  seasonal  variability  if  currents  were  to  become  weaker 
in  the  winter.  This  does  not  appear  likely,  though,  since 
Pouloin  and  Niiler  [1989]  found  relatively  small  eddy  kinetic 
energies,  representing  all  seasons,  of  about  110  cm2  s-2  in 
the  area  south  of  34*N.  This  is  consistent  with  a  real  south¬ 
ward  decrease  in  current  variability.  It  therefore  seems  likely 
that  the  trend  of  lower  energy  toward  the  south  is  primarily 
geographical,  rather  than  seasonal. 

All  told,  the  Eulerian  statistics  confirm  the  highly  vari¬ 
able  nature  of  the  California  Current  system.  The  core  jet 
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locations  differed  considerably  between  1987  and  1988.  Fur¬ 
ther,  in  most  cases,  the  standard  deviation  of  currents  is 
greater  in  magnitude  than  the  mean.  This,  of  course,  implies 
that  an  Eulerian  mean  California  Current  will  look  very  lit¬ 
tle  like  an  instantaneous  realization.  It  is  thus  interesting  to 
ask  whether  the  California  Current  has  a  sharp,  meandering 
high-velocity  core  throughout  the  entire  north-south  extent 
of  the  domain.  This  is  clearly  the  case  north  of  36*N,  but 
the  weaker  means  and  fluctuations  toward  the  south  could 
be  due  either  to  the  core  becoming  diffuse  and  weak,  or  to 
meandering  over  a  broad  region  of  otherwise  weak  variabil¬ 
ity.  In  order  to  gain  some  further  insight  about  this  point,  we 
found  the  strongest  individual  velocity  measurement  in  each 
0.5*  box  (Figure  5).  Large  velocities  are  found  at  all  lati¬ 
tudes,  so  the  meandering  jet  model  seems  to  apply  through¬ 
out.  There  is,  however,  a  tendency  for  the  peak  velocities 
to  be  about  20%  weaker  south  of  33*  N  than  north  of  that 
latitude.  This  north-south  decrease  in  core  velocities  coin¬ 
cides  roughly  with  the  location  where  both  the  mean  and 
fluctuating  currents  (Figure  4)  become  noticeably  weaker, 
and  Geosat  estimates  of  sea  surface  height  variability  de¬ 
crease  sharply  [Flament  et  at.,  1989].  Thus,  we  envision  the 
California  Current  as  a  meandering  but  coherent  jet  over 
the  domain  of  study,  with  a  core  velocity  maximum  which 
decreases  south  of  about  33*N.  The  kinematic  implication 
of  this  is  that  the  core  spreads  out  into  a  somewhat  broader, 
less  energetic  flow  as  the  current  proceeds  southward. 
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Fig.  5.  Velocity  corresponding  to  the  maximum  speed  observed 
in  each  0.5*  box  (total  data  set). 


4.  &>DY  DimisivinES 

4-1.  Single  Particle  Statistics 

Single  particle  dispersion  statistics  were  computed  in  the 
same  manner  used  by  Davis  [19856].  These  measure  the 
tendency  for  an  individual  water  parcel  to  deviate  from  the 
ensemble  average  trajectory.  Specifically,  the  eddy  diffusiv- 
ity  is  defined  as 

*""*(£)  =  ±  (la) 

where 

0  =  rn'(£,  t)rm'(£,  t),  (16) 

where  n  and  m  can  be  either  x  or  y,  r"(£,  ()  is  the 
location  of  a  drifter  in  the  n  direction  (since  the  time  of  its 
launch  at  fsO  at  position  £),  and  r"  is  the  difference 
in  the  n  direction  between  an  individual  drifter’s  location 
and  the  mean  drift,  rn  =  r"  —  (1  /N)  r".  The  average  is 
over  all  N  drifters  which  appear  within  a  particular  box. 
Since  the  memory  of  drifters  tends  to  be  short  (2-6  days;  see 
section  5),  all  drifters  were  treated  as  newly  deployed  after 
10  days  if  they  remained  in  the  box.  That  is  to  say  that, 
for  the  purpose  of  calculating  pnm,  a  single  drifter  could 
yield  several  trajectories,  each  starting  10  days  after  the 
other  if  the  drifter  was  in  the  box  at  its  new  starting  time. 
This  is  simply  a  means  of  retrieving  more  information  from 
the  available  data,  e.g.,  Poulain  and  Niiler  [1989].  For  all 
diffusivity  calculations,  the  coordinate  system  is  rotated  28* 
counterclockwise  so  that  t  is  perpendicular  to  the  general 
trend  of  the  California  coastline  between  34*30'N  and  39*N 
(Figure  1).  This  had  the  effect  of  decreasing  the  off-diagonal 
elements  of  the  dispersion  tensor  K**  in  some  cases,  while 
not  increasing  it  substantially  for  locations  farther  from  the 
coast. 

The  single  point  dispersion  calculations  were  carried  out 
for  three  boxes  having  particularly  good  data  coverage  (Fig¬ 
ure  1).  The  first  box  includes  the  filaments  within  which 
most  drifters  were  deployed.  The  second  box  is  farther  off¬ 
shore  and  includes  the  region  in  which  most  1988  drifters 
cease  going  offshore.  The  third  box  lies  in  the  southern  re¬ 
gion  where  the  drifters  tend  to  be  well  dispersed.  The  results 
for  the  merged  (1987  and  1988)  measurements  are  shown  in 
Figure  6.  In  all  cases,  the  xx  (onshore-offshore)  diffusivity 
is  the  largest  component  during  the  first  4  days.  For  the  two 
northern  boxes,  however,  the  diffusivity  decreases  by  day  5 
and  then  passes  through  zero.  This  behavior  may  reflect  the 
influence  of  long-lasting  flow  structures  (e.g.,  the  filament). 
For  the  southernmost  box,  no  such  persistent  flow  structures 
were  apparent  and  the  xx  diffusivity  had  more  of  a  ten¬ 
dency  to  increase  monotonically  with  time,  as  expected  in  a 
turbulent  field  [e.g.,  Davis,  19856]. 

In  all  cases,  the  initial  alongshore  yy  diffusivity  tended  to 
increase  with  time,  and  not  pass  through  zero.  This  suggests 
that  the  flow  structures  complicating  the  xx  dispersion 
act  mainly  in  the  onshore-offshore  direction,  as  might  be 
expected  for  meanders  oriented  roughly  perpendicular  to  the 
coast.  The  off-diagonal  diffusivity  was  small  for  t  <  10  days 
but  was  of  the  same  magnitude  as  the  alongshore  component 
of  diffusivity. 

The  results  of  the  single  point  diffusivity  calculations  are 
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Fig.  6.  Single  point  diffusivities,  calculated  with  the  merged 
(1987  and  1988)  data  sets.  Solid  line,  box  1;  dotted  line,  box 
2;  dashed  line,  box  3.  All  diffusivity  values  are  divided  by 
10TcmJ  s-1.  («)  Onshore-offshore  (xx)  compone'  ..  (i)  Along¬ 
shore  (yy)  component,  (e)  Off-diagonal  (xy)  component. 


summarized  (Table  1)  in  terms  of  the  average  values  over 
the  first  10  days.  Except  for  the  relatively  sparse  data  base 
in  1987  (±40%),  the  95%  confidence  intervals  on  the  10- 
day-averaged  diffusivities  are  typically  ±20%.  The  meaning 
of  traditional  confidence  intervals  in  this  situation  is  a  bit 
problematic,  especially  in  box  1,  since  the  drifters  are  not 
obviously  statistically  independent.  This  follows  because 
they  are  all  launched  into  an  energetic  flow  feature  which  was 
long-lived  relative  to  drifter  transit  times  of  5-10  days  [e.g., 
Kotro  et  al.,  this  issue].  Because  the  drifters  are  all  sampling 
the  same  feature  with  a  short  (1-2  months)  time  window, 
it  seems  likely  that  individual  drifters  are  not  statistically 
independent  and  that  the  means  and  diffusivities  do  not 
truly  reflect  an  ideal  “ensemble.” 

It  is  worthwhile  to  ask  whether  the  diffusivity  can  be  pa¬ 
rameterized  simply  in  terms  of  the  eddy  kinetic  energy  of 
the  system, 

where  v  is  the  vector  average  velocity  within  a  box.  Note 
that  since  this  measure  accounts  for  deviations  from  both 
the  spatial  and  temporal  mean,  it  will  tend  to  be  larger  than 
the  eddy  kinetic  energy  computed  at  a  single  site  from,  say, 
a  current  meter  mooring  (since  the  single  point  variance  is 
computed  relative  to  a  temporal  mean  only).  In  order  to 
reduce  scatter,  a  “scalar  diffusivity” 

A'  =  KIZ  +  Ayy 

averaged  over  10  days  lag  is  compared  with  E'K  (Fig¬ 
ure  7).  The  open  circles  correspond  to  values  calculated 
in  the  larger  boxes,  while  the  dots  correspond  to  estimates 
made  in  smaller  boxes  (dashed  boxes  in  Figure  1).  The  re¬ 
sult  suggests  a  dependence  of  K  on  E'h-  of  the  form 

K  =  b(E'h-?'\ 

where  b  =  42  km  (by  least  squares  fit).  This  result  sug¬ 
gests  [AVauss  and  Boning,  1987]  that  there  is  a  constant 
eddy  length  scale,  but  that  the  Lagrangian  time  scale  is 
not  constant,  a  finding  consistent  with  our  time  scale  cal¬ 
culations  (see  below  and  Table  1).  Our  isotropic  mixing 
length  of  42  km  can  be  compared  to  h’rauss  and  Boning’s 
[1987]  range  of  31-39  km  for  the  North  Atlantic  Ocean  (us¬ 
ing  their  definition).  Such  comparable  values  suggest  that 
our  measurements  correspond  to  an  open-ocean  regime,  not 
influenced  by  boundary  effects  in  a  gross  sense. 

The  averaged  single  point  diffusivity  is  always  initially 
largest  (for  short  time  lags)  in  the  onshore-offshore  direc¬ 
tion  (the  exception  being  in  1987  when  the  current  meanders 
were  less  pronounced).  To  interpret  this  result,  A’nm  can 
be  thought  of  as  the  diffusivity  to  be  used  in  a  steady  state 
diffusive-advective  model.  The  net  flow  is  southward  (or  to¬ 
ward  negative  y)  in  general  (Figure  4),  while  the  onshore- 
offshore  component  shows  substantial  variability  but  little 
net  trend.  Thus  advection  moves  the  drifters  fairly  consis¬ 
tently  toward  the  south,  while  the  “diffusion”  (parameter¬ 
izing  the  various  onshore  and  offshore  flow  patterns)  leads 
mainly  to  a  spreading  in  the  onshore-offshore  direction. 

4-2.  Taylor’s  Hypothesis 

If  the  flow  field  is  homogeneous  and  stationary  within 
a  box,  then  by  Taylor’s  hypothesis  the  eddy  diffusivity  is 
identically  equal  to 
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TABLE  1.  Eddy  Diffusivitie*  and  Independence  Time  Scale* 


10-Day  Average 

Single  Point  Kinetic  Independence 

Diffusivitie*,  Energy,  Time  Scale*,* 

cm3  *-*  x  10~T  cm2  i~3  _ day* 


Year 

Box 

nm»« 

K„ 

K „ 

Mean 

Eddy 

n 

n 

1987 

1 

43 

3.25 

1.88 

1.05 

114 

341 

2.1 

1.8 

1988 

1 

83 

10.04 

2.66 

2.32 

381 

811 

4.2(2.8) 

0.9 

2 

175 

5.45 

4.95 

-1.66 

56 

719 

2. 5(1.0) 

1.7(1. 7) 

3 

329 

5.61 

1.77 

-1.21 

100 

228 

5 .8(3.3) 

2.7 

Merged 

1 

126 

8.62 

2.43 

1.87 

227 

714 

3.3(2.8) 

0.9(0.8) 

2 

181 

5.34 

4.95 

-1.57 

54 

718 

2.5(1. 0) 

1.8(1. 6) 

3 

329 

5.61 

1.77 

-1.21 

100 

228 

5.8(3.3) 

2.7 

*  Maximum  number  of  drifter  days  in  the  calculation. 

t  From  (3)  and  (6).  Values  in  parentheses  are  time  scales  integrated  to  50  days-  Such  values  are 
only  given  if  the  result  stabilized  for  the  last  20  days.  Values  not  in  parentheses  are  for  integrations 
terminated  at  the  first  zero  crossing  of  the  autocorrelation  function. 


where  Lnm  is  the  lagged  correlation  function  of  the  La¬ 
grangian  n  and  m  components  of  flow  (t»n,ttm),  and 
un‘um‘  is  the  zero  lag  covariance.  The  Lagrangian  correla¬ 
tion  function  Lnm  was  computed  as 

N 

Lnm(t')  =  [tf<r"(0)<rm(t')]  £  u"'(0)um'(t'),  (3) 

where  the  standard  deviations  <rn  are  recomputed  sep¬ 
arately  for  each  time  lag  because  the  Lagrangian  mean 
changes  with  lag,  a  finding  consistent  with  Davis  [19856]. 
The  sum  is  over  all  N  drifter  trajectories  originating  in 
the  box.  Here  we  calculate  the  right-hand  side  of  (2)  for 
each  of  the  three  boxes.  Note  that  using  this  form  of  the 
autocorrelation  function,  which  is  necessary  to  account  for 
the  inhomogeneity  of  the  flow,  means  that  the  results  of  (1) 
and  (2)  need  not  be  identical.  Thus  the  extent  to  which 
this  estimate  agrees  with  the  Knm  estimated  using  (1)  is 
a  measure  of  the  homogeneity  of  the  flow. 

Comparisons  of  the  two  estimates  of  diffusivity  for  boxes  1 
and  3  are  shown  in  Figures  8.  Generally,  the  two  estimates 
agree  well  for  the  first  day  aft^r  drifter  launch,  but  substan¬ 
tial  differences  occur  within  d_ys.  The  divergence  can  be 
understood  through  examination  of  velocity  covariance  as 
a  function  of  time  since  launch.  The  cases  where  the  two 
estimates  diverge  quickly  correspond  to  case  where  the  co- 
variances  (e.g.,  kinetic  energy)  change  rapidly  as  a  function 
of  time  since  launch.  That  is,  divergence  of  the  estimates 
corresponds  to  cases  of  strongly  inhomogeneous  flow  statis¬ 
tics  over  scales  comparable  to,  or  smaller  than,  the  box  size. 
The  cases  where  the  two  estimates  diverge  slowly,  on  the 
other  hand,  correspond  to  situations  where  the  flow  statis¬ 
tics  remain  relatively  constant  with  time  after  launch.  It  is 
worth  noting,  too,  that  in  box  3,  Taylor’s  hypothesis  works 
fairly  well,  and  this  is  a  region  of  relatively  uniform  Ederi-n 
statistical  properties  (Figure  4). 

Taken  together,  the  tests  of  Taylor’s  hypothesis  lead  to 
the  conclusion  that  the  flow  in  the  region  sampled  off  Cali¬ 
fornia  is  rather  inhomogeneous  on  scales  as  small  as  200  km. 
This  should  perhaps  not  be  surprising  given  the  presence  of 


a  coast  and  the  apparent  north-south  variation  in  Eulerian 
statistics. 


4-3.  Particle  Pair  Dispersion 


The  rate  at  which  two  particles  separate  is  the  appropri¬ 
ate  eddy  diffusivity  to  consider  when  examining  the  time- 
dependent  behavior  of  a  system  [Davis,  19856].  This  particle 
pair  diffusivity  is  defined  by 


«nm 


1  £ 

4  dtl 


(4a) 


where 

*/nm  =  i£yVB'.  (46) 


Fig.  7.  “Scalar  diffusivity”  as  a  function  of  eddy  kinetic  energy. 
Circles  indicate  large  box  estimates  (solid  boxes  in  Figure  I). 
Dot*  indicate  small  box  estimates  (dashed  boxes  in  Figure  1). 
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Fig.  8.  Comparison*  of  single  point  diffusivitie*  computed  directly  (solid  line)  and  using  Taylor's  hypothesis  (dotted 
line).  All  diffiiaivity  values  are  divided  by  10ran3  s~* .  (s)  On»hore-off»hore  (rx)  component,  box  1.  (4)  Along¬ 
shore  (yy)  component,  box  1.  (c)  Onshore-oflshore  (xx)  component,  box  3.  (4)  Alongshore  (yy)  component, 
box  3. 


The  n  component  of  separation  between  two  drifters  is 
and  tn  is  the  deviation  of  a  given  drifter  pair  from 
that  of  the  mead  at  that  time  step.  The  average  is  over  all 
drifter  pairs,  with  no  geographical  binning.  Following  Davit 
[1985b],  the  calculation  of  vnm  was  binned  by  the  initial 
separation  distance  between  drifters.  No  attempt  was  made 
to  carry  out  separate  calculations  for  different  geographical 
areas  due  to  the  need  to  accumulate  large  numbers  of  de¬ 
grees  of  freedom,  although  nearly  all  pairs  started  in  either 
box  1  or  2  (Figure  1).  As  in  the  single  particle  dispersion 
calculations,  all  drifters  were  treated  as  newly  launched  af¬ 
ter  10  days.  The  bias  corresponded  to  initial  separations  of 
less  than  16  km,  between  16  and  32  km,  and  between  32  and 
64  km. 

The  results  of  the  calculations  for  the  merged  1987-1989 
data  set,  displayed  as  a  function  of  pair  separation  (which  in¬ 
creases  with  time  for  the  results  given)  are  shown  in  Figure  9. 
The  initial  separation  bin  can  be  identified  by  the  location 
where  a  given  curve  passes  through  or  »  0.  As  in  the  case  of 
•ingle  particle  diffusivitles,  the  onshore-offshore  zx  disper¬ 
sion  is  the  largest,  yielding  values  substantially  larger  (for 
the  same  separation)  than  those  reported  by  Davit  [1985b] 


for  the  nearshore  region  near  38*N.  The  alongshore  diffusiv- 
ities  are  far  smaller  than  the  onshore-offshore  component, 
apparently  as  a  result  of  the  alongshore  flow  component 
tending  always  to  have  an  equatorward  component,  as  would 
be  found  in  a  meandering  jet.  The  off-diagonal  components 
of  pair  dispersion  also  tend  to  be  small  and  show  little  net 
pattern.  In  all  cases,  the  dispersion  values  tend  to  decrease 
toward  zero,  corresponding  to  a  state  where  particle  pair 
separations  become  relatively  constant,  at  least  for  a  time. 
The  pair  diffusivities  tend  to  be  the  same  magnitude  as  the 
single  particle  dispersions,  although  pair  yy  values  tend  to 
be  somewhat  smaller  than  single  particle  yy  dispersions. 

5.  Dynamics 
5.1.  Integral  Time  Scalet 

As  Davit  [1985b]  points  out,cthe  autocorrelation  function 
of  velocity  can  be  a  useful  diagnostic  of  the  flow.  The  Eule- 
rian  independence  time  scale 

Tim  J~  A'Wt  (5) 


14,702 


Bunk  et  al.:  I'feAR-SuwxcE  Flow  in  CbiAsriAL  Transition  Zone 


g  xx  Pocr  Ovffuftlv\.ty 


yy  Pair  Diffuse vity 


xy  Pair  Oiffusivity 


12-1  £""(< 


(«) 


is  a  measure  of  how  quickly  the  velocity  ch  anges  for  a  particle 
being  advected  through  the  flow  field.  A  highly  nonlinear 
flow  field  dominated  by  semipermanent  flow  features  will 
likely  have  Tl  <  Te,  whereas  a  flow  field  dominated  by 
linear  dynamics  (flow  speeds  much  less  than  typical  “wave” 
speeds)  will  generally  have  Tl  —  Te  [Davit,  1985b]. 

In  order  to  make  such  a  comparison  for  the  coastal  tran¬ 
sition  cone  data  base,  the  Eulerian  time  scale  was  estimated 
from  the  current  meter  data  obtained  at  depths  of  145-200  m 
(data  described  by  Stabeno  and  Smith  [1987];  see  also  Fig¬ 
ure  1).  The  moored  current  meters  give  a  typical  analogous 
(one-sided)  value  of  12  days.  This  value  should  be  used  with 
caution,  however,  since  the  measurements  are  from  a  depth 
range  which  is  much  deeper  than  the  near-surface  bound¬ 
ary  layer  (see  section  3).  Indeed,  the  current  meter  velocity 
fluctuations  are  a  good  deal  weaker  than  the  surface  cur¬ 
rents  (typically  a  factor  of  4-5).  One  might  expect  that 
near-surface  velocity  measurements  would  yield  somewhat 
shorter  time  scales  due  to  the  presence  of  directly  wind- 
driven  flows,  but  estimates  of  momentum  balances  given  be¬ 
low  suggest  that  this  component  of  flow  may  not  be  large. 
Thus  the  145-200  m  current  meter  estimates  may  well  be 
reasonable  for  comparison. 

The  Lagrangian  time  scales  were  estimated  by  calculat¬ 
ing,  for  all  drifters  originating  in  a  given  box,  the  Lagrangian 
autocorrelation  function  by  (3).  All  calculations  were  car¬ 
ried  out  in  the  same  rotated  coordinate  system  as  was 
used  for  the  diffusivity  calculations.  Given  the  short  time 
scales  involved,  all  drifters  were  treated  as  newly  launched 
10  days  after  their  first  appearance,  leading  to  a  greater 
number  of  degrees  of  freedom.  In  many  cases,  the  integral 
(equation  (6))  did  not  converge  even  after  integrating  over 
50  days,  so  the  integration  was  instead  truncated  at  the  time 
of  the  first  zero  crossing  of  the  autocorrelation  function. 

The  results  of  the  calculations  (Table  1)  yield  Lagrangian 
time  scales  shorter  than  the  Eulerian  time  scale  (12  days), 
and  often  much  shorter.  Thus  the  hypothesis  of  fundamen¬ 
tally  nonlinear  dynamics  appears  to  be  validated,  at  least  in 
the  north  where  current  meter  data  are  available.  Notice, 
however,  that  the  idea  of  nonlinear  dynamics  does  not  mean 
that  the  nonlinear  terms  in  the  momentum  equation  need 
to  be  large.  The  result  means  only  that  such  terms  must  be 
important  relative  to  the  other  ageostrophic  terms. 


5.2.  Momentum  Balance* 

It  is  possible  to  use  the  drifter  data  to  estimate  the  terms 
in  the  upper  ocean  momentum  budget.  Specifically,  we  write 
the  horizontal  momentum  equations  as 

AS  ,  ,  .  1_  ldf  .  . 

3j  +  yi».  +  -v,---  =  o  (7) 


Fig.  9.  Pair  diffusivities  divided  by  10rcm*  »~* .  Initial  sep¬ 
aration  bins  can  be  distinguished  by  where  the  curves  have 
o  =  0.  (<)  Onshore-offshore  (xx)  component,  (i)  Alongshore 
(w)  component,  (c)  Off-diagonal  (xy)  component. 


(where  A "  is  the  autocorrelation  function  of  the  n  com¬ 
ponent  of  velocity)  is  a  measure  of  the  typical  time  scale 
of  velocity  changes  for  an  observer  located  at  a  fixed  point. 
The  analogous  Lagrangian  time  scale 


where  small-scale  horizontal  turbulent  mixing  has  been  ne¬ 
glected.  The  total  derivative  is 


dt 


d  .  d  .  d  ^  3 

—  4*  u 4- 1? —  4-  it; — , 

at  dx  dp  a*' 


(») 


where  /  is  the  Coriolis  parameter,  k  is  the  vertical  unit 
vector,  p  is  the  density  of  water,  and  f  is  the  vertical 
component  of  turbulent  stress.  If  a  slablike  mixed  layer  is 
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HINO  STRESS  AT  46013 


DAYS  IN  1988 

Fig.  XO.  Wind  strew  at  buoy  46013.  North  ia  upward.  For  reference,  July  1  it  day  183.  Horizontal  line*  represent 
the  times  of  hydrographic  surveys. 


assumed  [e.g.,  Davit  et  al.,  1981],  then  the  vertical  gradient 
of  turbulent  stress  in  the  mixed  layer  can  be  stated  as 

~  =  A"1*.,  (9) 

where  A  is  the  mixed  layer  depth  and  ft  is  the  wind  stress 
applied  at  the  ocean’s  surface. 

Using  the  drifter  data,  it  is  straightforward  to  estimate  the 
Coriolis  force  term  in  (7).  Further,  the  total  derivative  of 
velocity  is  also  straightforward  to  estimate,  except  that  the 
drifter  can  not  account  for  the  vertical  derivative  term  in  (8). 
The  problem  then  is  to  estimate  the  pressure  gradient  and 
the  mixed  layer  depth.  The  approach  taken  was  to  consider 
only  those  drifters  in  or  near  the  hydrographic  sampling  grid 
(Figure  1)  during  the  time  when  a  particular  hydrographic 
survey  was  made,  typically  a  period  of  6-8  days.  Then, 
during  this  time,  the  first  two  terms  in  (7)  were  regressed 
against  the  surface  wind  stress,  so  that  the  regression  co¬ 
efficient  a  could  be  identified  with  (pA)"1.  There  were 
four  hydrographic  surveys  during  which  there  were  sufficient 
drifters  to  carry  out  this  calculation  (surveys  3-6).  The  first 
three  occurred  during  fairly  strong  winds,  and  the  last  dur¬ 
ing  weak  winds  (Figure  10).  In  only  one  case  was  a  reason¬ 
able  value  of  A  obtained:  map  3,  where  A  =  26  m  but 
the  x  correlation  is  only  0.17.  In  all  other  cases,  either  a 
negative  or  large  O(>80  m)  positive  value  was  obtained  for 
A,  in  contrast  to  observed  mixed  layer  depths  of  0(20  m). 
Finally,  the  pressure  gradient  was  calculated  as  the  resid¬ 
ual  of  (7)  and  all  results  were  averaged  into  0.5*  by  0.5* 
(latitude  by  longitude)  Inns  to  clarify  presentation. 

Typical  results  of  a  calculation  are  shown  in  Figure  11, 
representing  survey  3  on  days  187-194  (July  6-13),  1988. 
The  Coriolis  (arrow)  and  acceleration  (short  stub)  terms  are 
compared  ia  Figure  11a,  while  the  stress  gradient  and  Cori¬ 
olis  terms  are  compared  in  Figure  lli.  Clearly,  the  Coriolis 
terra  dominates,  suggesting  a  nearly  geostrophic  balance. 


The  drifter  estimate  of  pressure  gradients  (dashed  arrow)  is 
compared  to  the  pressure  gradient  based  on  0/499  dbar  dy¬ 
namic  height  in  Figure  11c.  There  is  a  qualitative  agreement 
between  the  two  fields,  but  it  is  not  very  good.  The  momen¬ 
tum  balance  results  from  the  three  other  surveys  yielded 
comparable  results. 

The  poor  results  of  the  regression  approach  should  not  be 
surprising.  For  a  velocity  within  the  filament  of  50  cm  s_I , 
a  Coriolis  term  of  0(5  x  10~s  cm  s-2)  would  be  antici¬ 
pated.  For  a  typical  wind  stress  of  1  dyn  cm”2  and  mixed 
layer  depth  of  20  m,  a  typical  Ekman  velocity  would  be 
about  5  cm  s-1  and  the  stress  gradient  term  in  (7)  would 
be  0(0.5  x  10_I  cm  s-2),  small  enough  to  play  only  a 
secondary  effect  relative  to  the  dominating  geostrophic  bal¬ 
ance,  and  making  the  regression  approach  difficult.  The 
poor  agreement  between  the  two  pressure  gradient  esti¬ 
mates  could  have  two  causes:  (1)  errors  in  computing  the 
momentum  balance  from  drifters  (including  the  fact  that 
the  drifters  are  not  perfect  Lagrangian  instruments),  and 
(2)  499  dbar  not  being  the  correct  level  of  no  motion  for 
the  flow  field.  The  moored  current  results  of  Stabeno  and 
Smith  [1987]  suggest  that  a  correct  level  of  no  motion  in  this 
area  would  lie  at  a  depth  of  1000  dbar  or  greater  and  that 
errors  associated  with  the  494-dbar  reference  level  would  be 
0.0003-0.0012  cm2  s-1,  too  small  to  account  for  all  of  the 
discrepancy  in  Figure  11  A,  but  still  limiting  our  ability  to 
assess  the  quality  of  the  momentum  balance  estimates  from 
comparisons. 

Nonetheless,  some  conclusions  can  be  drawn  from  the  mo¬ 
mentum  balance  calculation.  First,  the  lo west-order  mo¬ 
mentum  balance  must  be  geostrophic.  This  finding  is  consis¬ 
tent  with  the  validity  of  quasigeostrophic  dynamics.  At  the 
next  order,  Ekman  transport  and  other  ageostrophic  effects 
probably  enter.  It  is  not  possible  to  say  directly  whether 
the  total  time  derivative  is  dominated  by  local  or  advective 
changes.  Results  of  the  time  scale  calculations  (section  5.1) 
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Fig.  11.  Remit*  of  momentum  balance  calculation*  for  day*  188- 
198, 1988.  («)  CorioB*  (with  arrowhead)  and  acceleration  (with¬ 
out  arrowhead)  term*.  Location*  are  given  by  dot*.  (1)  Cori- 
oli*  (with  arrowhead)  and  (tree*  gradient  (without  arrowhead) 
term*,  (c)  Comparison  pressure  gradient*  computed  from  drifter* 
(dashed  arrow*)  and  from  0/499  dbar  dynamic  height  (solid  line*). 
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Fig.  1 1 .  (continued) 


suggest,  however,  that  the  nonlinear  terms  must  be  impor¬ 
tant. 

6.  Relation  of  Near-Surface  Flow  to  Fronts 

Very  often,  filaments  have  fronts  along  both  their  north¬ 
ern  and  southern  sides.  Further,  there  is  some  evidence  that 
convergence  and  subduction  could  be  occurring,  at  least  at 
the  southern  front  [Flamcnt  et  al.,  1985].  It  thus  seems 
reasonable  to  attempt  to  use  the  drifter  data  to  look  for 
evidence  of  such  convergence. 

The  approach  used  is  straightforward.  Frontal  locations 
are  determined  subjectively  by  using  NOAA  AVHRR  satel¬ 
lite  thermal  imagery  for  1988  only.  These  are  overplotted 
with  drifter  positions,  so  that  the  shortest  distance  of  a  given 
drifter  from  both  the  southern  side  (ij)  and  northern  ride 
(()  fronts  can  be  determined.  Note  that  this  approach  does 
not  depend  upon  absolute  surface  temperatures,  but  rather 
only  on  the  distance  from  fronts.  Such  features  should  not, 
over  short  time  scales  (days),  depend  greatly  on  surface  heat 
fluxes.  Most  drifters  were  deployed  within  the  cold  waters  at 
the  nearshore  end  of  the  Point  Arena  filament.  Those  that 
were  initially  not  between  the  fronts  are  dropped  from  the 
calculation.  Histograms  of  the  distance  of  drifters  from  the 
front  versus  time  since  launch  were  then  computed  (sum¬ 
marised  as  means  and  standard  deviations  in  Figure  12). 
As  the  drifters  proceed  offshore  in  the  filament,  they  draw 
closer  to  the  south  side  front,  i.e.,  if  decreases.  The  av¬ 
erage  distance  to  the  north  side  front  (  remains  roughly 
constant  at  about  15  km.  This  is  consistent  with  the  fila¬ 
ment  narrowing  with  distance  offshore  as  is  typically  seen 
in  satellite  imagery.  These  results  are  also  consistent  with 
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Fig.  12.  («)  Average  distance  of  drifters  from  the  south  side  front 
(rj)  versus  time  since  launch.  Vertical  bars  represent  ± 1  standard 
deviation.  (1)  Average  and  standard  deviation  of  drifter  distance 
from  the  north  side  front  «).  (e)  Average  velocity  out  the 

filament  axis  versus  time  from  launch. 


convergence  on  the  sonth  side  of  the  filament,  but  they  do 
not  mean  that  down  welling  (subduction)  is  occurring.  For 
example,  the  offshore  flow  in  tne  filament  could  simply  be 
accelerating  to  compensate  for  the  horizontal  constriction. 

The  acceleration  hypothesis  was  tested  by  computing  the 
average  velocity  as  a  function  of  time  for  all  of  the  drifters 


used  above.  The  result  (Figure  12)  shows  that  initially  the 
drifters  did  accelerate  but  that  they  slowed  down  by  day  5 
from  launch.  Thus  even  though  the  filament  becomes  nar¬ 
rower,  the  flow  along  its  axis  decelerates,  implying  down- 
welling.  Most  likely,  this  downwelling  occurs  at  or  near  the 
southern  side  front,  since  the  drifters  tend  to  converge  on  it 
(Figure  12).  It  is  straightforward  to  compute  a  volume  bal¬ 
ance  for  the  flow  over  the  upper  15  m  of  the  water  column 
(assuming  no  flow  through  the  filament  edges)  to  arrive  at 
a  downwelling  rate  of  12  m  d-1  in  the  section  within  the 
filament  3  to  5  days  after  deployment  (100-170  km  from  de¬ 
ployment).  This  rate  of  sinking  is  roughly  consistent  with 
the  20  m  d“*  estimate  based  on  biological  and  chemical 
tracers  [Kadko  et  al.,  this  issue]. 

The  drifter-based  sinking  velocity  estimate  is  not  very 
precise.  It  involves  many  assumptions,  including  that  no 
flow  crosses  the  edge  of  a  filament.  Further,  errors  are  likely 
to  be  introduced  by  the  unavailability  of  satellite  data  dur¬ 
ing  cloudy  periods,  which  lasted  up  to  5  days.  Thus,  the 
results  need  to  be  taken  with  caution.  This  is  especially  so 
because  Swenson  et  al.  (1990)  concluded  that  downwelling 
takes  place  on  the  northern  (not  southern)  edge  of  the  fila¬ 
ment.  While  their  result  does  not  directly  call  into  question 
our  estimate  of  vertical  velocity,  it  does  serve  as  a  caution 
about  the  idea  that  downwelling  tends  to  occur  preferen¬ 
tially  along  the  southern  side  of  the  filament. 

7.  Conclusions 

The  CTZ  drifter  results  strongly  suggest  that  the  Califor¬ 
nia  Current  between  30*  and  40*N  has  a  meandering,  jetlike 
current  core.  Some  of  the  structure  in  the  current  (e.g.,  the 
eastward  flow  near  36*N)  appears  to  be  reproducible  from 
1987  to  1988,  although  the  fairly  featureless  mean  patterns 
south  of  34*N  suggest  that  the  core’s  position  varies  fairly 
randomly  in  that  area.  The  magnitude  of  the  core  velocity, 
the  mean  velocity  and  the  eddy  kinetic  energy  all  tend  to 
decrease  toward  the  south,  especially  south  of  about  33*N. 

Although  the  lowest-order  momentum  balance  for  the 
near-surface  currents  is  geostrophic,  the  governing,  higher- 
order  dynamics  are  more  complex.  Specifically,  the  disparity 
between  Eulerian  and  Lagrangian  current  time  scales  sug¬ 
gests  that  nonlinear  effects  are  quite  important.  This  finding 
is  consistent  with  the  interpretation  of  the  California  Cur¬ 
rent  system  as  being  hydrodynamically  unstable,  e.g.,  Ikeda 
and  Emery  [1984]  and  Haidvogel  et  al.  [this  issue]. 

Accompanying  the  nonlinearity  of  the  flow  is  a  tendency 
for  rapid  particle  dispersion.  The  values  for  both  the  single 
point  and  pair  dispersions  are  substantially  larger  than  those 
estimated  by  Davis  [19856]  for  an  area  closer  to  shore  near 
38.5*N.  Davis  [19856]  also  noted  that  dispersion  tended  to  be 
isotropic  by  50  km  from  shore,  while  alongshore  dispersion 
dominated  close  to  the  coast.  Our  results,  which  represent 
the  region  greater  than  50  km  from  shore,  show  a  strong 
tendency  for  dominance  by  the  onshore-offshore  component 
of  dispersion.  This  is  the  opposite  sense  of  anisotropy  to  that 
found  by  Davis  [19856].  Presumably  this  difference  is  related 
to  the  coastal  inhibition  of  onshore  flow  over  the  shelf,  and 
to  the  dominance  of  meandering  jet  patterns  offshore. 

The  flow  field  offshore  is  sufficiently  inhomogeneous  as  to 
limit  Taylor’s  hypothesis  to  only  order  of  magnitude  valid¬ 
ity.  This  inhomogeneity  is  also  reflected  in  the  substantial 
spatial  variations  in  kinetic  energy  and  eddy  dispersion  co- 
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efficients.  This  inhomogeneity  is  roughly  consistent  with 
sea  level  variability  estimated  with  Geos  at  [Flament  et  al., 
1989]. 

It  now  seems  very  likely  that  subduction  occurs  in  fila¬ 
ments  off  northern  California,  based  on  chemical  and  hy¬ 
drographic  tracers  [Kadko  et  al.,  this  issue].  The  results  of 
the  CTZ  drifter  study  suggest  that  much  of  this  occurs  more 
than  100  km  offshore,  and  at  an  area-averaged  rate  within 
the  filament  of  0(10  m  d-1).  Apparently  the  cold,  dense 
water  of  the  filament  core  eventually  has  to  sink  below  the 
warmer  ambient  water.  By  the  time  the  meander  causing 
the  filament  reverses  back  toward  shore,  the  cold  surface  sig¬ 
nature  is  thus  greatly  reduced.  This  would  account  for  why 
satellite  temperature  and  color  imagery  only  tend  to  “see” 
flow  proceeding  ofTshore. 
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